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fibers 



O 

o 

(N 

>>: 



00 



Shaoping Bian, Dirk Robinson, and Mark G. 
Department of Physics & Astronomy, Washington State University, Pullman, Washington 99164-2814 

(Dated: February 2, 2008) 

We report on what we believe is tlie first demonstration of an optically activated cantilever due 
to photomechanical effects in a dye-doped polymer optical fiber. The fiber is observed to bend 
when light is launched off-axis. The displacement angle monotonically increases as a function of 
the distance between the illumination point and the fiber axis, and is consistent with differential 
light-induced length changes. The photothermal and photo-reorientation mechanisms, each with its 
own distinct response time, are proposed to explain the observed time dependence. The measured 
degree of bending is consistent with a model that we have proposed which includes coupling between 
photoisomerization and heating. Most importantly, we have discovered that at high light intensity, 
a cooperative release of stress results in cis-to-trans isomerization that yields a large and abrupt 
length change. 
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I. INTRODUCTION 



While the intensity-dependent phase shift of light in 
a material is often studied, change in the length in a 
mode of propagation in the material is not usually con- 
sidered; rather the intensity dependent refractive index 
and absorption mechanisms are assumed to dominate. 
However, there is a long history of observations of the 
photomechanical effect. Indeed, in 1966, Merian re- 
ported that fabrics doped with azo-dyes contract under 
light exposure. Eisenbach found that a pre-stretched 
chromophoric poly(ethyl acrylate) network shrinks along 
the direction of aUgned polymer chains when exposed 
to light. More recently, reversible shape changes in 
solids in the range from 10% to 400% and in films up 
to 20% that are optically induced by photoisomeriza- 
tion of monodomain nematic elastomers were shown by 
FinkelmannJJ and Li et al 3, respectively. Camacho- 
Lopez and coworkers showed that small eleastomer sam- 
ples that float on water will "swim" away from regions 
exposed to light. @ The strong coupling between orien- 
tational order and material strain was described by de 
Gennes in the early 1970s. 

In this paper, we report the observation that light, 
which is propagating down a fiber, can be used to induce 
it to bend, thus functioning as an optically activated can- 
tilever. Many mechanisms of intensity dependent length 
change have already been discussed. 0, ll| Using the 
photothermal mechanism,'^ several devices have been 
demonstrated including a tunable optical filter, jlO] an 
all-optical vibration suppressor jl^, 0| and a mesoscale 
version of such a device that exhibits both mechanical 
and optical multistabilitv. [T^ In this paper, we show that 
differential expansion can be used to make an all-optical 
cantilever in a dye-doped polymer optical fiber, and we 
study its mechanisms through modelling and experiment. 
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We begin by reviewing some of the more rudimentary 
observations. This leads us to the development of our 
coupled theory of the photomechanical effect - which orig- 
inates from photothermal heating - and photoisomeriza- 
tion. This section is followed by a set of experiments that 
focus on measuring the parameters of the theory, such as 
the intensity dependence of the time constants and mag- 
nitude of bending as a function of pump intensity. We 
find that at low intensities, the theory does not behave 
as expected, because the length decreases rather than in- 
creasing with pump intensity. However, above a thresh- 
old intensity, the length change becomes abruptly posi- 
tive and large - which is reminiscent of a process in which 
internal stress builds until a threshold is reached, after 
which the stress is released. So, we propose a hypothe- 
sis that collective interactions between the molecules are 
responsible. We find that our results are consistent with 
this model; but, we should keep in mind that the model 
is not meant to be a precise theory. Rather, its purpose 
is to deduce the underlying physics of the mechanisms of 
the response. 



II. LASER INDUCED BENDING OF A FIBER 

When a laser beam is launched off-axis into a short 
fiber made of poly(methyl methacrylate) (PMMA) that 
is doped with Disperse Red 1 azo-dye (DRl), the fiber is 
observed to bend. Fig.^a) shows a micrograph of a such 
fiber when no laser beam is launched into it. The fiber 
end and sides are aligned to a reticle grid as shown with 
the three arrows prior to illumination. Fig. ^b) shows 
a micrograph of the same part of the fiber when a laser 
beam is launched in from the right, and the arrow on the 
right of Fig. ^h) indicates the launch position. Relative 
to the superimposed grid, the fiber's length is observed 
to increase by about half of a division and bending re- 
sults in the fiber end moving downwards by almost a 
full division. This cantilever is found to move back and 
forth reproducibly between the two states when the light 
is repeatedly turned on and off; so the phenomenon is 
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FIG. 1: A photograph of an MPU (a) without iUumination 
and (b) with iUumination. The three arrows mark reference 
positions for identifying changes of the fiber position. The 
white arrow indicates the position of the excitation beam. 



reversible. 

The bending is caused by differential expansion be- 
tween the illuminated and dark portions of the fiber. The 
illuminated side expands more than the darker side, and 
the degree of the bending depends on the distance of 
the pump beam from the fiber axis. Fig. |21 shows the 
measured degree of bending as a function of the pump 
beam's displacement from the fiber axis. As expected, 
the bending angle increases linearly as a function of the 
axial beam position when the pump beam is fully con- 
tained within the fiber. When the pump is near the fiber 
edge, some of the light propagates outside the fiber, lead- 
ing to a decrease in the amount of bending. 



III. THEORY 

In this section, we develop a photomechanical theory of 
heating and photo-reorientation mechanisms. This cal- 
culation is not meant to be complete in that it does not, 
for example, treat the full three-dimensional problem. 
Rather, the goal is to provide a semi-quantitative frame- 
work that is applied to understanding the observations. 
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FIG. 2: The measured degree of fiber bending as a function of 
the position of the pump laser, measured from the fiber axis. 
The inset shows the fiber end and pump beam, approximately 
to scale. 



To this end, several simplifying approximations are made. 
These approximations can be understood in terms of the 
following description. 

Consider a polymer optical fiber with an embedded 
isotropic distribution of one-dimensional chromophores 
(such as DRl) in the trans state that undergo polarized- 
laser induced isomerization to a cis state that is bent, and 
therefore smaller. The orientational mobility of the cis 
molecules will be much larger than the trans ones. The 
trans molecules, oriented along the polarization direction 
of a light beam, will have the highest probability of op- 
tical absorption. As such, molecules that are oriented 
parallel to the laser's polarization axis will be converted 
to the cis isomer, which subsequently will reorient due 
to the higher degree of mobility and then decay back to 
the trans state. The net result is that molecules oriented 
parallel to the laser's polarization will be depleted, and 
converted to an anisotropic distribution, with the net re- 
sult of a "hole" in the orientational distribution function. 
Consequently, more molecules will be oriented perpendic- 
ular to the laser polarization than parallel to it. 

We assume that there are two mechanisms of light- 
molecule interactions that lead to photomechanical 
effects. In the photo-reorientation process, photo- 
isomerization results in a net re-distribution of molecules 
in a plane that is perpendicular to the light's plane of 
polarization, which includes a component that is perpen- 
dicular to the fiber axis and another component along the 
fiber axis. This leads to a stress along the fiber axis re- 
sulting in a commensurate change in length and a change 
in the fiber diameter perpendicular to the beam's po- 
larization, producing an elliptical fiber cross-section. In 
the photothcrmal process, the absorbed energy is turned 
into heat, which results in thermal expansion. Since the 
photo-isomerization process causes an orientational hole, 
the amount of energy absorbed by the fiber drops over 
time, so the photothcrmal process efficiency drops. In 
this way, the two mechanisms are coupled. 
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As a zeroth-order approximation, we assume that 
photo-isomerization is independent of the temperature 
for the range of temperature changes expected for photo- 
heating. In this case, the photo-reorientation mechanism 
can be calculated independently of photothermal heating. 
The resulting orientational distribution of chromophores 
due to the photoisomerization process can then be used 
as a source term in the photothermal calculation, thereby 
coupling the two mechanisms. Secondly, we will ignore 
the effect of the change in fiber diameter and focus only 
on the change in fiber length. As such, we will assume 
that the length change and ellipticity change are decou- 
pled; and, that the ellipticity change does not result in 
bending. Indeed, a small part of the cantilever effect that 
we observe could be attributed to changes in fiber ellip- 
ticity; but, if it is proportional to the length change, such 
an effect could be taken into account as a small correction 
factor in proportion to the length change. 



A. Photo-Reorientation Mechanisms 

A full calculation of the mechanisms would require us 
to consider a general orientational distribution function, 
which even for the simple case of the photo-reorientation 
mechanism alone yields complex results that require nu- 
merical integration. We will therefore simplify the prob- 
lem by effectively using a two-state model; that is, we 
assume that a molecule can be oriented only in two or- 
thogonal orientations: perpendicular or parallel to the 
light beam's polarization. We will also ignore the cis 
isomer population and assume that it is small at any in- 
stant of time relative to the population of the resulting 
reoriented trans molecules. 

We note that a three-state model could be a better ap- 
proximation because the molecules that are reoriented 
away from the polarization direction are equally dis- 
tributed in the perpendicular plane. Since on average, 
we can view such a distribution as half of the molecules 
being aligned along the fiber axis - resulting in a length 
change, and the other half being oriented along the fiber 
diameter - resulting in a diameter change, a three-level 
model might be more appropriate. Such an effect can be 
simply taken into account by attributing a length change 
to only half the population that is not aligned along the 
light's polarization direction. Thus, a fit of the two-level 
theory to experiment should yield a length change pa- 
rameter that is double the value that would be given by 
a three-state model; but should yield the correct dynam- 
ics. So, we proceed in our theoretical development with 
the understanding that the parameters of our theory will 
correctly account for the dynamics of our observations; 
but, that a more rigorous calculation would be required 
to relate these parameters to real physical properties of 
the polymer and dopant chromophores. 

The approximations we use are as follows: 

1. The trans molecule interacts with light only if it is 
oriented along the polarization of the light beam. 



2. The trans molecule, after interacting with a photon, 
leads to a trans molecule oriented perpendicular to 
the photon's polarization. 

3. When the material is not illuminated, an cntropic 
process causes the system to relax into a steady 
state equilibrium with equal populations of both 
orientations. 

We define the following variables: 

• N - The fraction of molecules oriented along the po- 
larization of the light beam. 1 — iV is thus the frac- 
tion of molecules perpendicular to the polarization; 
and, = 1/2 is the non-illuminated population. 

• / - The intensity of the light. 

• ^ - The probability per unit of intensity per unit of 
time that a trans molecule will absorb light and be 
converted to the perpendicular orientation. 

• /3 - The entropic decay rate of an angular hole. 
The population, A^, can thus expressed as 

dN 



dt 



-S,IN + P{l-2N), 



(1) 



where (1 — 2A^) is the difference between the perpen- 
dicular and parallel population fraction. Eq. 1^ can be 
integrated to yield. 



A^ : 



2/3 + g/cxp[- (2/? + g/) t] 
2 (2/3 + e/) 



(2) 



where we have evaluated the integration constant by de- 
manding that N{t = 0) = 1/2. Note that we can also 
express this result in a way that better shows the form 
of the exponential growth. 



A^^l 
2 



1 - 



(2/3 



(l-exp[- (2/3 + e/)t]) 



(3) 



There are several points about this result that need 
to be mentioned. First, at infinite time, the equilibrium 
population is: 



eq 



2/3 + 



(4) 



and the conversion rate of reorientation is intensity de- 
pendent and of the form. 



6(1) = 2/3 + e/. 



(5) 



The limiting forms of Eq. are reasonable. When the 
entropic decay rate is much larger than the rate of light- 
induced reorientation (/3 ^I), then N^q = 1/2 as ex- 
pected. If the decay rate /3 is small compared with the 
reorientational conversion rate then all of the parallel 
population is converted and N^q = 0. 
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When the Hght source is turned off, if the initial pop- 
ulation is N^q and the decay rate is given by Eq. ^ 
with / = 0, this yields a population N oriented along the 
pump beam polarization of 



N 



exp(-2/3t), 



(6) 



2(3 + ^1 



where we have used Eq. Note that Eq. (|SJ) yields 

N{t ^oo) = 1/2. 



B. Photothermal Heating Mechanisms 

In this section, we calculate the temperature as a func- 
tion of time under the assumption that heat is trans- 
ferred from the light beam to the sample through optical 
absorption of only the parallel population. If the light 
beam is turned off, we assume that the temperature fol- 
lows Newton's law of cooling. We define the following 
parameters: 

• To - The temperature of the substance surrounding 
the fiber. 

• 7 - The cooling rate of the fiber. 

• a - The temperature increase per trans molecule 
parallel to the beam's polarization per unit of in- 
tensity per unit of time. Clearly, a depends on the 
concentration of dopants. 

The heating rate of a fiber under uniform illumination 
is then given by. 



dT 



-7(r - To) + aNI, 



(7) 



where the first term is the Newton cooling rate and the 
second term represents energy deposited through optical 
absorption. When the light is turned off, the cooling 
process is described by: 



(8) 



Clearly, the constant 7 depends on the fiber geometry 
and its specific heat. 

a). Light-on illumination process 

We begin by considering the heating process. Substi- 
tuting Eq. (0) in Eq. O and using Eq. lO, we get: 



dT , , al 



l-^^[l~expM.)] 



(9) 



Integrating Eq. 0, we get 
a/5/ 



T-To - Aexp(-7i)- 



al [8 - 2/3) 



2 (7 - 5) 



exp(— (5t), 
(10) 



where the first term on the righthand side is the homoge- 
nous solution with integration constant A and the rest of 
the expression is the inhomogeneous term. To find A we 
demand that at i = the fiber is in thermal equilibrium 
with its surroundings so that T = Tq. This yields. 



A = - 



af3I 



al {S - 2/3) a/(7 - 2/3) 



(57 2 (7 -(5) S 27(7 -(5) 

Substituting Eq. (HTJ into Eq. yields: 



(11) 



T-To 



al 



2(7 - S) 
5-213 



7-2/5 
7 



[1 - exp(-7i)] (12) 



[1 - cxp(-(5t)] 



b). Light-ofF ^ decay process 

Next, we solve for the cooling process, give by Eq. (jS)). 
This yields. 



T-To = Bexp(-7t), 



(13) 



where B is an integration constant. The typical experi- 
mental sequence consists of turning on the light source, 
waiting until the fiber reaches equilibrium (i.e. t — > 00), 
then turning off the light source. As such, we choose the 
temperature difference T—Tq at i = as the temperature 
difference from Eq. \i2\ at infinite time, or 



T-To = 



aj3 



I = B, 



so the final cooling result is 



^-^0 = (^)/exp(-7t). 



(14) 



(15) 



C. Photomechanical Response 

We are now ready to formulate the photomechanical re- 
sponse. First, we consider the photo-reorientation mech- 
anism. When a molecule reorients away from the laser's 
polarization axis, the strain along that axis should de- 
crease and the strain perpendicular to it should increase. 
So, the fractional length change, AL/L (where L is fiber 
length and AL is its change caused by the photomechan- 
ical effect) should be proportional to the difference be- 
tween the fraction of molecules oriented away from the 
polarization and those parallel to it. 



AL 
17 



= b[{l - N) - N] =6(1 -27V), 



(16) 



where h is what we call the isomer-mechanical constant, 
which depends on the properties of the polymer and 
dopant molecules. Note that the length changes if the 
dopant molecules deviate from their isotropic distribu- 
tion of iV = 1/2. 
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FIG. 3: A schematic diagram of the relationship between the 
length change and degree of bending. 



In the photothermal mechanism, the length changes 
due to thermal expansion. For a coefficient of thermal 
expansion at, and assuming additivity between the two 
mechanisms, the total length change is given by 



^ = h{l-2N) 



oitiT - To) , 



(17) 



Next, we need to determine the relationship between 
the fractional length change and bending angle of the 
fiber. As stated above, if the temperature is higher on 
one side of the fiber, that side will expand, leading to 
bending away from the hot side. We assume that half 
the fiber cross-section is uniformly illuminated while the 
other half is uniformly dark, so that issues such as con- 
tinuous temperature gradients are ignored. Fig. |21 illus- 
trates a bent fiber along with the parameters we use in 
the following calculation. If the fiber is of length L be- 
fore bending, we assume that after bending, the length 
of the fiber along the non-illuminated side remains L. If 
the length increase of the illuminated side of the fiber is 
AL, we see from Fig. O that 



R + d 



L ^ LAL AL 



where e = L/d is the length-to-diameter ratio of the can- 
tilever, and R is the radius of curvature as shown in Fig- 
ure |31 

We can calculate the dynamics of photomechanical 
bending of a fiber by substituting AL/L in Eq. H18|) into 
Eq. (|17() using the appropriate expression for the number 
fraction iV and the temperature difference T — Tq from 
the previous section. 

It is important to note the effects that we have ne- 
glected in the above derivation. First, we have assumed 
that half the fiber cross-section is uniformly illuminated 
and the other half is dark. Furthermore, we assumed 
that the temperature increase is proportional to the lo- 
cal intensity, that there are no field gradients, and the 
the cooling rate depends on the fiber being in contact 
with a heat batch (i.e. air). 

Clearly, the fiber is not illuminated uniformly since 
the beam has a gaussian profile. In addition, parts of 



the beam undoubtedly reflect from the fiber/air inter- 
face making the intensity profile even more complicated. 
Secondly, a temperature gradient will form between the 
illuminated region and the dark side, which complicates 
the temperature profile and can also cause the beam to 
refract into the dark region. These effects together lead 
our calculations to overestimate the amount of bending 
and to underestimate the time constant. Furthermore, 
since our approximation assumes that the heated side 
follows its natural length and that the dark side follows 
suit, this also overestimates the bending angle. However, 
our beam is focused so that a region that is smaller than 
half the cross-section is illuminated, which should result 
in a larger temperature gradient, thus leading to an un- 
derestimate of photo-bending. Furthermore, we assume 
that the cylindrical fiber has a rectangular cross-section. 
Since the magnitude and time response of the amount of 
bending that we observe for a polymer fiber agree well 
with our calculation, we conclude that for the level of 
accuracy needed to understand the dynamics of the re- 
sponse, our simplified model is sufficient. Even if a better 
calculation were available, there are too many unknown 
parameters for such a model to be of any utility. As such, 
we proceed with our simplified model of photothermal 
heating. 

a). Light-on illumination process 

We begin by considering the pump illumination pro- 
cess. Substituting Equations 10) and H12I) into H17|l . and 
the result into Eq. (|18|l . the change in angle due to the 
photoisomerization and photoheating mechanisms, 6^ , is 



eatal (7 - 2/3) 
27 (7 - 6) 



- exp(- 
atal 



'It)] 



2(7-5) 



[1 - exp(-(5t)] (19) 



Note that we have combined terms with the same time 
constant. Finally, we substitute Eq. |SJ| into Eq. (|19|l to 
get the result in terms of the material parameters, 



eatal (7 - 2/3) 
27 (7 - 2/3 -e/) 



[1 - exp(-7i)] 



Thermal Rise - Fast 



2(3- 
(1- 



exp I 



atal 



^ "2 (7 - 2/3 - il) 
-(2/3 + e/)i]) . 



(20) 



Orientation Rise - Slow 

Each process has associated with it a characteristic time 
scale. As we will see below, photoisomerization is the 
slower process, so we label that term the slow process 
and the heating process term the fast process. 

b). Light-ofF^ decay process 

Similarly, for the decay process when the pump illumi- 
nation is off, we substitute Eq. ijHJl and H15|) into Eq. H17|) 
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FIG. 4: A schematic diagram of the experiment used to mea- 
sure bending of the MPU in response to light excitation. VA: 
variable attenuator, M: mirror, L: lens, and CL: cylindrical 
lens. 



and the result into Eq. I|18(l . This yields the angle, 6"^, as 
a function of time in a dark fiber, 



eat a/3 J 



exp(~7t) - 



eb^I 

2/3 + a 



exp(-2/3t). (21) 



Thermal Decay Orientation Decay 

Note that the initial value of Eq. H21|) (at i = 0) is equal 
to that of Eq. idni) at t = oo. 



IV. EXPERIMENT 

To systematically study the photomechanical effects, 
an experimental setup schematically shown in Fig. ^ is 
adopted. The fiber used in this experiment is made of 
PMMA that is uniformly doped with 1% (by weight) 
DRl dye. The diameter of the fiber is 600 /im and the 
length is 2.6 mm, yielding an aspect ratio of e = 4.33. 
The fiber's absorption coefficient is a = 0.3 mm~^ at a 
wavelength of A = 633 nm, at which the excitation laser 
operates. We call such fibers mesoscale photomechani- 
cal units (MPUs). Details of the fabrication process for 
making an MPU can be found in the literature. 0, 0, 0) 

The MPU cantilever is polished on both ends, and 
one end is attached to a glass substrate with transpar- 
ent cyanoacrylate adhesive. A small piece of a metalized 
glass cover slip is bonded to the other end and functions 
as a lightweight mirror. A beam from a He-Ne laser, 
called the pump beam, is coupled into the fiber though 
the substrate. The pump laser beam is focused by a lens 
and the MPU is located just beyond the focal point, so 
the spot size of the laser on the input end face of the 
fiber is estimated to be ^ 150 /im. An electronic shut- 
ter is used to control or modulate the light. A second 
laser beam, which we call the probe beam, is incident on 
the MPU's mirror and the reflections from the mirror is 




FIG. 5: A schematic diagram showing how a length change 
or bending of an MPU affects the measured power at the 
detector. 



directed onto an opaque plate with a narrow slit behind 
which a photodetector is placed. A cylindrical lens is in- 
serted between the MPU mirror and detector (see Fig. 
to expand the probe beam on the opaque plate thereby 
increasing the spatial resolution of the slit/detector pair. 
The pump beam is set off axis near the edge of the fiber 
for maximum deflection. Also, the pump laser beam is 
offset horizontally from the fiber axis so that the induced 
MPU bending occurs in the plane of the table (i.e. in 
the plane of the page in Fig. ^ and therefore the move- 
ment of the probe laser beam is perpendicular to the slit 
on the opaque plate. Both pump and probe lasers are 
He-Ne lasers operating at A = 633 nm. 

Fig. |S1 shows how the change in length and bending 
affects the beam. Clearly, the offset of the beam due to 
a length change is independent of the distance between 
the MPU's mirror and the detector while bending leads 
to larger deflections when the distance is made larger. 
The change of the detected power due to bending thus 
depends on this distance and the intensity profile of the 
laser beam. The greatest sensitivity is achieved when 
the slit is placed at the point of highest intensity gra- 
dient in the beam. Furthermore, the dynamic range of 
deflections that can be measured depends on the beam 
diameter. Dynamic range is therefore inversely related to 
the sensitivity. In our experiments, the intensity distri- 
bution of the probe laser beam is Gaussian. We choose 
the beam size, focal length of the cylindrical lens, power 
range of the pump laser, and the slit position such that 
only the left half (expanded by the cylindrical lens) of 
the Gaussian beam translates across the slit for the full 
range of laser-induced MPU deflection. The detector is 
placed far enough away from the cantilever to increase 
the sensitivity of detecting bending while minimizing the 
effect of length change. The output of the detector is fed 
to a digital oscilloscope in which the experimental data 
can be saved. In our setup, we have verified that the de- 
tected power change due to the MPU's length change is 
much smaller than that caused by bending of the MPU. 
Therefore, the power change measured by the detector 
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FIG. 8: The deflection angle of the probe laser (circles) as a 
function of time after the shutter is opened. The gray curve 
is a fit to Eq. lHJ. 
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FIG. 7: Degree of bending of an MPU cantilever as a function 
of time with pump laser repeatedly turned on and off. 
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function of time after the shutter is closed. The gray curve is 
a fit to Eq. 



reflects mainly the effects of bending. 

The apparatus is calibrated to determine the de- 
gree of angular deflection per unit of power change of 
the probe beam. To do so, the MPU assembly (sub- 
strate/MPU/mirror) is placed on a rotation stage with 
a rotation angle resolution of 0.017°. The sample is ro- 
tated in small, well calibrated angular increments and the 
power is measured in units of detector voltage. Fig. El 
shows the data and a polynomial fit. This calibration 
polynomial is used to convert the measured probe beam 
power into a bending angle of the cantilever. Note that 
above a deflection angle of 0.1°, the displacement of the 
beam spot at the detector plane is larger than the beam 
diameter, and therefore beyond the detectable range of 
the slit/detector combination. As such, our experiment 
is limited to deflection angles of less than 0.1°. 



V. RESULTS AND DISCUSSION 

The temporal response of the cantilever is obtained by 
measuring the power at the detector as a function of time 
after the shutter is opened or closed. Fig.EJshows the out- 
put of the detector when the shutter is opened and closed 
repeatedly, which illustrates that the fiber moves repro- 
ducibly back and forth. Fig.|Slshows a typical plot of the 
deflection angle (converted from the measured power at 
the detector by using the calibration polynomial obtained 
from the fit shown in Fig. (S)) as a function of time after 
the shutter is opened for a 12.7 W/ cm? pump intensity. 
Here we have used an estimated average pump beam di- 
ameter of 200 /iTO inside the MPU. We call the process 
while the shutter is open a rise process. Fig. |^ shows a 
continuation of Fig. |S1 when the shutter is closed, which 
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FIG. 10: The equilibrium deflection angle of the MPU as a FIG. 11: Time constants for both the (a) fast and (b) slow 
function of pump power process as a function of pump power obtained after the pump 

is turned on. Each point in the plot is determined from data 

fits similar to Fig.|S| 



we call a decay process. The characteristics of these re- 
sponses are as follows. They can be well fit to Eq. H20() 
and H21(l respectively, with two different time constants, 
as shown by the gray curves in Figs. |H1 and El (all fit pa- 
rameters are shown as insets). This indicates that there 
are two mechanisms that are responsible for the response 
as modelled by our theory. Furthermore, the two expo- 
nential components have opposite signs (there are a few 
exceptions in the rise data, which we will discuss below). 
According to Eq. H20() and (|21|l , it can be concluded that 
at low intensity. 



6 < 0. 



(22) 



This implies, according to Eq. (|16|l . that the bending 
component caused by the photo-reorientation mechanism 
is negative, or, photoisomerization yields a decrease in 
the fiber length. 

The deflection angle is measured as a function of time 
after the shutter is opened; and, once equilibrium is 
reached, after the shutter is closed, at several pump in- 
tensities in the range of 3.2 W/cm'^ to 16.9 W/cm^. Equa- 
tions 1)21) |l and (|21|l are used to fit these data. As such, 
we obtain the power dependence of the physical param- 
eters in our theoretical model, such as the equilibrium 
steady-state bending angle (Fig.^J; the rise and decay 
time constants for the fast and slow mechanisms fFies.lTTl 
and ll2() : and the amplitude of the rise and decay process 
for the fast and slow mechanisms fFig. 1131 and I14|l . In 
Fig. El for intensities above YlW/cnri? , the deflection an- 
gles were measured by observing the reflected spot on a 
screen and determining the angle trigonometrically. Any 
Intensity above \lW/cm? corresponds to a deflection an- 
gle greater than 0.1°, which is beyond the range of the 
calibrated slit/detector pair. 
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FIG. 12: Time constants for both the (a) fast and (b) slow 
process as a function of pump power obtained after the pump 
pulse is turned off. Each point in the plot is determined from 
data fits similar to Fig. |n| 



A. Time Constants 

According to our theory, the rate, 7, for heating and 
cooling are the same and independent of the intensity for 
the rise and decay processes (see Eq. (|20f) and i'Zll ). In 
contrast, the rates of photo- reorientation and its relax- 
ation are different in the rise and decay process, because 
the former is power dependent (i.e. due to the term 
2/3 -I- ^/ in Eq. H2UII ) and the latter is power indepen- 
dent (see Eq. idU). Therefore, the hypothesis that the 
two time constants are different for the photoisomeriza- 
tion processes is relatively simple to test with intensity- 
dependent studies of the slow mechanism. 

By comparing the data in Figs. ^2 and El we find 
that within experimental uncertainty, only a weak in- 
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FIG. 13: The fast components of the amplitude of the deflec- 
tion angle as a function pump power obtained from (a) rise, 
and (b) decay data. The open circles indicate the ratio of 
these two components. 



FIG. 14: The slow component of the amplitude of the deflec- 
tion angle as a function pump power obtained from (a) rise, 
and (b) decay data. 



tensity dependence is observed. Under some conditions, 
response times on the order of seconds in the photo- 
reorientation phenomenon have been reported in the 
literature. [11 J5, 16, 17, 18] However, since the time con- 
stant depends on the material composition and intensity, 
it is not a straightforward matter to determine the time 
constant that one would expect in our experiments. 

To determine which time constant corresponds to 
which mechanism, we calculate the thermal time con- 
stant of the MPU using data from our experiments and 
PMMA's thermal properties from the literature. 19] The 
calculation takes into account the geometry and mass of 
the fiber and the material's thermal properties (such as 
the specific heat), and the result indicates a time con- 
stant of about a few hundred milliseconds for our experi- 
mental conditions, which is consistent with the fast time 
constants we obtained (shown in Fig. IllT a) and I12r a)). 
Therefore we infer that the fast process originates from 
photothermal heating, and the slow process is due to 
photo-reorientation. 

An average of the response times in each set of data 
in Fig. ^] and E| gives a mean fast response time of 
Ty — 0.113 (±0.005) s from the rise process and = 
0.109 (±0.003) s from the decay process; while for the 
slow process, we get rj — 3.36 (±0.30) s from rise data 
and rf = 2.83 (±0.15) s from decay data (these numbers. 



except tJ, are labelled in Fig. 1111 and ll2() . The fact that 
the fast rise and decay time constants are independent of 
each other and are equal within experimental uncertainty 
supports our hypothesis that the fast process must be of 
a photothermal origin. 

In the following section, we extract the various param- 
eters of our model by analyzing the dependence of the 
data on time and intensity. 



1. Heating — Fast Mechanism 

From the data, an average rise and decay rate of 7 = 
1/t/ = l/(0.11s) = 9.01 (±0.16) s-i is found. 



2. Photo-Reorientation — Slow Mechanism 

a). Decay process 

The decay rate due to the entropic process of relax- 
ation to the isotropic state is power independent, and is 
determined from Eq. to be 
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b). Rise process 

The rise time constant originates from the competition 
between intensity-induced excitation and entropic decay, 
so is power dependent. Using the relationship between 
the response time and the intensity from Eq. H20|l . the 
rise time constant is 



1 



2/3 + a 



(24) 



Fitting the data in Fig. Illf b) to this func- 
tion yields (3 = 0.13 (±0.04) s"! and C 
0.0044 {±0.0078) cm'^W-^s-\ The value of (3 obtained 
here from the rise data is 28% smaller than that ob- 
tained from the decay data. However, they agree with 
each other within experimental uncertainty. An average 
of/3 = 0.16s~^ will be used in following calculations. 
But, we found that scatter in the rise data is relatively 
larger, and of the fit, which is the sum of the squares 
of deviations of the theoretical curve from experimental 
points, is 0.95. Therefore, the values obtained from this 
fit, especially the value of is subject to a large uncer- 
tainty. This uncertainty can be further revealed in the 
following analysis. 



B. Amplitude 

1. Heating — Fast Mechanism 

a). Rise process 

In our experimental conditions, we have 7 — 2/3 >> 
and from Eq. H20|) the bending angle of the fast rise 
process at t — > oo is 



negligible is not surprising because the data in Fig. lisr b) 
do not show any evidence of saturation, as predicted by 
Eq. (EEJ. Because the value of the fit is 8 13 x 10"^ 
it is much better than the fit in Section lV A 2b : therefore 
providing a more accurate value of ^ than we obtained 
in Section IV A2b . When ^ is neglig ible, Eq. ^ is the 
same as Eq. H25II . and both amplitudes of rise and decay 
process follow a linear dependence on pump intensity. 
The ratio of these two amplitudes calculated from our 
experimental data is shown in Fig. IKSf b) (open circles), 
which demonstrates that the experimental data is con- 
sistent with our theory. 



2. Photo-Reorientation — Slow Mechanism 

a). Decay process 

The slow decay amplitude given by Eq. (|21|l is 



9^(t = 0)= '^^^ . 



(27) 



Fitting Eq. H27|) to the experimental data in Fig. I14r a). 
we get = -1.9 (±0.3) x lO-^cm^W-^-^deg, and, 
again we find that ^ is negligible. Here we have used 
an average entropic relaxation rate of /3 = 0.16s~^, as 
obtained above. The value is 1.58 x 10^^, so the fit is 
good. Note that we kept the term ^/ in the denominator 
of Eq. H27|) in our fitting routine so that we could obtain 
the parameter ^ and compare it with the corresponding 
one we obtained in Section IV B lb . 

b). Rise process 

Since 7 — 2/3 >> £,1, according to Eq. H20|l . the ampli- 
tude of the photo-reorientation process is 



9}{t^ 00) 



eatal 



(25) 



Fitting the power dependence of the fast compo- 
nent of the bending angle shown in Fig. I13r a') to this 
equation gives the parameter ata = 1.94 (±0.03) x 
10^^ cm?W deg . Here we have used the length-to- 
diamctcr ratio e = 4.33 for our cantilever. 

b). Decay process 

According to Eq. (|21|l . the fast component of the decay 
of the bending angle at t = is 



-"^-7(2/3 + 



(26) 



Fitting the power dependence of the fast amplitude 
component shown in Fig. I13r b) to this equation gives 
the parameter ata = 2.02 (±0.01) x 10~^ cm^H^~^s~^def/ 
and a negligible value of ( < 10^^^ cm^W~^s~^). The 
value of the parameter a^a is in good agreement with 
the one obtained from the rise process. The fact that ^ is 



2/3 + e/ 



atal 



2 (7 - 2/3) 



(28) 



Since b < 0, 0^ should be negative. The experimental 
data, as shown in Fig.^Jb), are negative at low power, 
but positive at high power (i.e. the three open circles 
are of the same sign as the bending angle of the heating 
process). A hypothesis that explains the sign reversal 
will be presented below. If we fit the negative data of 
Fig.Etb) only to Eq. igHl, we obtain b£ = -1.2 (±0.1) x 
10~^ cm^W s^^ deg, which is about two-thirds of the 
value we obtained from the decay process. In this fit, the 
parameter ^, again, is found to be negligible; and so is 
the second term in brackets in Eq. H28|l in comparison 
with b. 

Our hypothesis for why the amplitude of the slow 
photoisomerization mechanism starts out negative and 
abruptly jumps to a large positive value above XZW jcrr? 
of pump intensity is as follows. We begin by arguing 
why a negative angle is observed. To illustrate the im- 
portant features of our physical picture, we consider only 
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FIG. 15: Each frame shows a microscopic view of those 
molecules that are affected most by the light (i.e. we ignore 
those molecules that are not along the light's polarization) 
and a macroscopic view of the fiber at a given applied field 
due to a light beam. The dashed lines show the original shape 
of the material or fiber while the solid line shows the instan- 
taneous shape. Note that Ei < E2 < E3. 



molecules that are aligned along the direction of the po- 
larization of the light beam since they are the ones that 
are excited with high probability. Furthermore, it is well 
known that the distribution of voids in a polymer vary 
over a large range of sizes. For this part of our argument, 
we only consider those voids that tightly fit around the 
molecule; and for convenience, represent them as ideal 
ellipsoids. Fig. Iisr a) shows an illustration of the ini- 
tial states that are under consideration. Both the fiber 
and microscopic view of the molecules are shown. While 
the axes of the molecules are randomly oriented so all 
orientations are equally represented, we only show those 
molecules in Fig.^lthat are along the light's polarization 
direction. Also, we show only those voids that are com- 
parable in size to the molecular size (since these are the 
ones that are responsible for the collective process that 
we describe below) even though larger voids are present. 
Note that we will discuss the role of larger voids at the 
end of this section. 

A negative angular deflection of the fiber is caused by 
a decrease in the length of the illuminated side where 
the polarization is perpendicular to the fiber axis, as 
shown in Fig. llSf b). The light therefore causes the 
molecules that are oriented perpendicular to the fiber 
(i.e. parallel to the light's polarization) to be excited to 
the cis state followed by randomization of the orienta- 





(d) 



FIG. 16: The process of trans-cis-trans molecular re- 
orientation in a void that is large enough to easily accom- 
modate the new orientation. 



tion of these smaller cis isomers. Subsequently, if any 
cis isomer decays back to a trans state that is oriented 
along the light's polarization, it causes those molecules 
to be re-excited. At higher intensity, a significant por- 
tion of the trans population becomes depleted as shown 
in Fig. ll5r c). Since the cis molecules sweep out a smaller 
volume than the trans ones, we would expect a decrease 
of void size around the cis molecules, accompanied by a 
decrease in the polymer's bulk volume, and therefore a 
decrease in the fiber's length. While we draw our fibers at 
high temperature, we have found that there is still a small 
amount of residual axial stress. 20j Shrinkage of a poly- 
mer in the direction of stress in response to illumination is 
commonly observed. To summarize, at low intensities 
(below liW/cm?), our model assumes that the process 
follows a sequence of events as shown in Fig. I15r a)-(c). 
This process is reversible. 

We present the hypothesis that the abrupt transition 
from negative to positive bend angle arises from a release 
of internal stress due to the population of cis molecules 
that are not otherwise conformationally able to relax 
to the trans state (as shown in Fig. llSf c)). While the 
smaller cis molecule may be able to rotate or tunnel into 
a state with its long axis perpendicular to the cylindrical 
void, it is difficult for such a cis molecule to decay into a 
trans state in this orientation. One can imagine that the 
cis molecule is being prevented from decaying into the 
trans state, leading to an internal stress in the polymer 
whose net direction is on the voids' narrow surface (left 
to right in Fig. llSf c)). We thus speculate that if enough 
molecules accumulate aligned in the same direction in 
the cis state, they cooperatively force the polymer to ex- 
pand along the narrow part of the void, allowing them 
all to collectively decay into the trans state as shown in 
Fig. Eld). So, the aggregate force due to the dyes yields 
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an increase in the fiber's length as shown by the solid 
lines in the diagram of the fiber in Fig. llSf d). Though 
the elongated polymer is in a higher energy state due to 
increased strain, the overall energy is lower due to the 
release of stress in the local environment of the interact- 
ing monomers and due to the presence of the light beam. 
This argument is analogous to ferroelectric domain align- 
ment where the energy of the system is lowered as the 
surface area of domain walls decreases. (However, in our 
model, residual stress remains, as described below.) 

If the collective alignment process is elastic in nature 
(i.e. the polymer remains stressed due to internal elastic 
forces and conformational changes because the viscosity 
is negligible on the time scales of the measurement), en- 
tropic processes cause fluctuations in alignment that can 
result in realignment, relieving the residual strain. The 
presence of the light, however, prevents enough chro- 
mophores from realigning collectively and lowering the 
energy, so the fiber remains in its elongated state. Only 
when the light beam is turned off can entropic processes 
lead to a collective reorientation accompanied by a de- 
crease of the fiber length. Since this process is observed 
to be fully reversible, we argue that the random state of 
orientation is the favored one, so the presence of the light 
beam plays an important role in driving the collective 
reorientational response that leads to a length increase, 
and without the presence of the light, the random state 
of orientation is of lowest energy. Furthermore, while we 
have focused on the illuminated side of the fiber, clearly, 
the dark side - which is stressed by the illuminated side 
- probably also play a role in making a straight fiber the 
most energetically favored state. 

The observation of a small decrease in length at low 
powers (Fig. llSf c)) and a sharp change to a positive de- 
flection angle at larger power (Fig. llSf dll is consistent 
with the picture that at low power, the polymer shrinks 
because there are not enough cis molecules to collectively 
overcome the internal stress barrier while at higher pow- 
ers, the barrier is exceeded. While illuminated, the length 
of the fiber is determined by a balance between the stress 
due to the polymer that acts to decrease the length and 
the collective force of alignment, which is driven by the 
light beam. When the source of illumination is turned 
off, entropic forces lead to a collective relaxation of chro- 
mophore alignment. 

For our hypothesis to be plausible, the concentration 
of DRl molecules must be large enough to allow cooper- 
ative interactions between them. For the 1% by weight 
concentration used in our studies, the mean distance be- 
tween nearest neighbors is about twice the length of the 
DRl molecule. As such, it is reasonable to expect that 
molecules near each other can mechanically interact us- 
ing the polymer as a mediator: As one void is deformed, 
the affects are felt by nearby voids as the polymer chains 
react by steric interactions. 

We must reconcile this observation with the literature 
on the measured refractive index change observed under 
similar conditions. The universal observation is that the 



refractive index change grows monotonically with time 
and as a function of the intensity. Furthermore, the 
measured birefringence shows that molecules are being 
converted continuously to trans isomers that are perpen- 
dicular to the polarization direction of the beam. This 
observation appears at first to be inconsistent with our 
observation of a discontinuity. We reconcile this appar- 
ent contradiction by arguing that the distribution of sites 
in polymers provides a means for the same microscopic 
processes to result in different observations for optical 
and mechanical behavior. (We note that given the ge- 
ometry of our experiment, it was not possible to measure 
the change in absorption spectrum of the illuminated side 
while operating the cantilever, so the optical properties 
are deduced from thin-film measurements.) 

The distribution of voids in a polymer is known to 
be large. As such, we would expect that some of the 
molecules are loosely held in place with voids whose 
sizes are approximately equal to the length of the long 
axis of the trans state while others are held tightly in 
place by the smaller ellipsoidal voids. The loosely-bound 
molecules can easily reorient to any direction in the trans- 
cis-trans sequence. Fig. llBf a) shows a trans molecule in 
such a void while Fig. Iltir b') shows how the cis molecule 
can freely rotate. So, we argue that these reoriented dye 
molecules as shown in Fig. I16f c) will yield a negligible 
stress, which can all be accommodated by small realign- 
ment of the chains near the molecule without a global 
change in the polymer's shape. So, reorientation can be 
observed optically, but leaves no mechanical trace. The 
tightly bound molecules, however, can only reorient if 
the local free volume changes shape drastically, which 
by virtue of the large conformational changes required of 
the polymer, can only occur if there are global changes 
in the polymer's shape. Such global changes require the 
collective action of many molecules as seen in Fig. lTST d'). 



VI. CONCLUSION 

We have demonstrated differential photomechanical ef- 
fects in dye-doped polymer fibers. The optically-induced 
bending of a fiber is governed by two mechanisms with 
different time response. The analysis suggests the fast 
one is attributed to a photothermal process whilst the 
slow one is shown to be related to polymer matrix de- 
formation caused by polarized-laser-induced ordering of 
the azo-dye molecules. More importantly, we find what 
appears to be a new phenomena that leads to photome- 
chanical effects: At high intensities, internal stress ap- 
pears to build until there are enough molecules to cause 
the polymer to reconfigure itself in a way that collectively 
releases the stress. 

Since our work has focused on the fundamental physi- 
cal basis of the all-optically actuated cantilever, we have 
intentionally studied systems with low dye concentra- 
tion because of the resulting simplifications in the mod- 
elling. A higher degree of bending is routinely observed in 
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beams, plates, and films |2 If using isomerization in li quid 
crystalline materials^ and shape-memory polymers. |22j 
There are potentially many important applications that 
may result from miniaturization of devices 23^ such as 
nano-positioning, optical switching. (24. .25] and beam- 
steering. The experimental techniques and theory pre- 
sented here can be used to study new materials, design 



devices, and predict their behavior. 
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